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* 

R.  E.  LeLevier 

The  RAHD  Corporation,  Santa  Monica,  California 

ABSTRACT 

A  theoretical  mydel  is  presented  for  calculating  the  heating 
and  hydrodynamic  expansion  of  the  F-layer  of  the  ionosphere  due  to 
an  incident  radiowrve .  We  find  that  the  greatest  effects  occur 

for  transmitter  frequencies  just  above  the  critical  frequency  of  the 
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F-layer.  Using  a  power  aperture  product  of  10  megawatt-meter  , 
which  corresponds  to  a  projected  experiment  by  ESSA  at  Boulder. 
Colorado,  we  find  a  maximum  increase  in  the  electron  temperature  at 
the  peak  of  the  F-iayer  of  about  thirty  percent  with  a  reduction  in 
the  electron  density  of  about  ten  percent.  The  time  scale  for 
achieving  the  density  changes  is  predicted  to  be  about  one-half  hour 
to  an  hour. 


♦Any  views  expressed  in  this  paper  are  those  of  the  author.  They 
should  not  be  interpreted  as  reflecting  the  views  of  The  RAND  Corporation 
or  the  official  opinion  or  policy  of  any  of  its  governmental  or  private 
research  sponsors.  Papers  are  reproduced  by  The  RAND  Corporation  as  a 
courtesy  to  members  of  its  staff. 

This  paper  was  prepared  for  publication  in  the  Journal  of  Geophysical 
Research . 


1  -  INTRODUCTION 


There  has  been  considerable  interest  in  the  problem  of  heating 
the  electrons  in  the  F-region  of  the  ionosphere  by  radiowave 
propagation.  Ginzburg  and  Gurevich  [1960]  and  Gurevich  [1967]  in  the 
Soviet  Union  have  published  extensively  on  this  problem  and  Farley 
[1963]  in  the  United  States  has  also  addressed  the  problem.  In  the 
work  done  to  date  steady-state  conditions  have  been  assumed  and,  in 
addition,  pressure  equilibrium  as  a  function  of  altitude  has  been 
assumed  for  the  heated  plasma.  These  approximations  are  reasonably 
valid  for  equatorial  latitudes  where  the  earth's  magnetic  field  lines  are 
more  or  less  horizontal.  With  horizontal  field  lines  the  heated 
plasma  can  equilibrate  horizontally  with  the  ambient  unheated  plasma. 

In  northern  latitudes  where  the  earth's  magnetic  fields  are  more 
nearly  vertical,  it  would  appear  that  the  assumption  of  pressure 
equilibrium  is  not  appropriate.  Indeed,  the  cross  field  drift  of 
the  plasma  due  to  lateral  pressure  gradients  is  so  slow  (less  than 
1  cm/sec)  that  the  heated  plasma  is  essentially  confined  to  expand 
up  and  down  along  the  field  lines.  We  are  interested  it;  this 
geometry  since  ESSA  plans  to  conduct  experiments  at  Boulder,  Colorado 
this  fall  in  modifying  the  normal  F-layer. 

In  the  following  section  a  theoretical  model  is  described  which 
approximates  the  hydrodynamic  motion  as  being  in  one  dimension  along 
the  field  lines.  Thermal  conduction  of  the  heated  electrons  is 
accounted  for  as  well  as  the  viscous  drag  imposed  on  the  ions  by  the 


2 


neutral  particles  of  the  ambient  atmosphere.  The  neutrals  are 
treated  as  a  heat  sink  at  the  ambient  temperature  of  the  ionosphere. 

The  partial  differential  equations  representing  the  conservation 

equations  were  numerically  integrated  on  a  high  speed  computer.  The 

re  nj  r-p  of  the  calculations  are  shown  in  Section  3  where  we  examine 

the  temperature  rise  and  the  rate  of  expansion  of  a  model  F-layer  as 

a  function  of  the  transmitter  frequency.  In  the  calculations  we 

4  2 

have  assumed  a  power-aperture  product  of  10  Megawatt  meter  , 
continuous  wave,  appropriate  to  the  GSSA  experiment.  We  find  that  a 
maximum  heating  rate  is  obtained  for  a  frequency  which  almost  gives 
a  null  for  the  dielectric  constant  of  the  ordinary  ray.  This  result 
differs  from  that  of  Gurevich  [1907")  who  finds  an  optimum  heating 
frequency  at  about  twice  the  critical  frequency.  This  difference 
may  be  traced  to  the  fact  that  Gurevich  maintains  a  constant  beam 
width  aa  the  transmitter  frequency  is  varied.  That  is,  Gurevich 
treats  the  problem  of  a  fixed  incident  flux  on  the  ionosphere.  This 
can  be  accomplished  by  changing  the  size  of  the  antenna  aperture  as 
a  function  of  frequency  or  by  varying  the  power  as  the  frequency 
changes.  We  have  considered  a  somewhat  more  practical  example  of 
constant  aperture  size  which  corresponds  to  the  planned  experiment. 
Therefore,  for  frequencies  well  above  the  critical  frequency,  we 
find  that  the  heating  rate  is  essentially  Independent  of  frequency. 
However,  as  one  reduces  the  frequency  and  approaches  the  null  in  the 
dielectric  constant,  an  enhancement  of  the  heating  rate  occurs.  As 
an  example  of  this  we  will  present  the  results  of  a  heating  program 
which  maintains  the  transmitter  frequency  to  within  IX  of  the  nulling 
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frequencies  as  a  function  of  time. 

We  will  begin  by  describing  the  tneoretical  model  and  then  will 
present  the  numerical  results. 
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2 .  THEORETICAL  MODEL 

We  imagine  that  an  antenna  is  directing  radiowave  energy  into 
the  ionosphere  and  parallel  to  the  earth's  magnetic  field  lines.  At 
the  latitude  of  Boulder,  Colorado  this  would  require  a  pointing  angle 
of  about  22°  from  the  zenith,  the  dip  of  the  magnetic  field  lines 
being  68°  at  that  latitude.  As  the  RF  energy  heats  the  plasma  of 
the  ionosphere,  the  plasma,  of  course,  tends  to  expand.  We  observe 
that  the  magnetic  pressure  is  far  greater  than  the  plasma  pressure 
and  the  cross-field  drift  of  the  plasma  due  to  the  pressure  gradient 
is  extremely  small.  Therefore,  the  expansion  of  the  plasma  will  be 
*mainly  upwards  and  downwards  along  the  field  lines.  Thus,  the 
familiar  approximation  of  pressure  equilibtium  [Farley,  1968;  Gurevich 
19671  is  not  appropriate  and  we  find  it  convenient  to  numerically 
integrate  the  conservation  equations,  thereby  investigating  the 
transient  approach  to  steady-state  conditions.  Since  we  do  not 
expect  unusually  large  departures  from  ambient  conditions,  we  will 
consider  a  perturbation  theory  approach  as  outlined  in  the  following. 
A.  Momentum  Equation.  We  consider  a  macroscopic  velocity  of  the 
plasma  and  write  the  momentum  equation  as 

0  Dt  +  +  6  "to"  +  0) 


where  the  plasms  density  p  is 
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p  *  m+N+  +  m  N  *  (m++  m  )N,  (2) 

m  and  m  are  the  ion  mass  and  electron  mass,  respectively,  N  is  the 
electron  number  density, 

6P  *  P-*  “  P+  +  P  *  Pa  (3) 

where  M+"  and  denote  quantities  referring  to  ions  and  electrons 
respectively,  and  the  subscript  "a"  denotes  ambient  values.  The  term 
°vioU  aPProxinUte8  the  "omentum  exchange  to  the  ambient  neutral  particles, 
vio  being  the  ion-neutral  collision  frequency.  The  last  term  in  Eq.  (1) 
is  due  to  the  external  force  of  gravity,  where  g  is  the  acceleration 
of  gravity,  *  is  the  dip  angle  of  the  earth's  field  and  60  is  the 
departure  of  the  plasma  density  from  ambient.  The  quantity  X  is  the 
distance  along  the  field  line  and  it  related  to  altitude  c  by 


X  sin*  *  *.  (4) 

Finally,  D/Dt  is  the  Lagrange  time  derivative  and  the  coordinate  X  is 
related  to  u  by, 


(3) 
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B.  Conservation  of  Mass.  Conservation  of  mass  in  the  plasma  is  given 
by 

0dX  -  0  dx  =  dM  (6) 

a 


where  M  is  the  usual  mass  variable  and  x  is  the  value  of  X  at  f«o. 

C.  Conservation  of  Energy.  We  treat  the  electron  gas  and  the  ion  gas 
separately  thereby  accounting  for  possible  differences  in  the  electron 
and  ion  temperatures.  If  we  define  the  specific  energy  of  the  two 
gases  by 


_+ 

IT— 


2  + 
m™ 


(7) 


where  0~  is  the  ion  (electron)  temperature  in  energy  units,  we  can 
write  the  partial  differential  equations  expressing  energy  conservation 

as 


m  N 


rq  .  i  i  i.  k  -j 

**  ax  *  & 


(8) 


where  K  is  the  coefficient  of  heat  conduction 


K  *  3He’/(m’v’) ,  (9) 

and  L  denotes  the  loss  of  energy  by  electrons  to  ions  and  neutrals. 

Electrons  lose  energy  to  both  nitrogen  molecules  and  atomic 
oxygen.  The  elec  tton  energy  loss  rate  to  nitrogen  molecules  arises 
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from  rotational  excitation  [Mentzoni  and  Row  1963].  However,  a  perhaps 
dominant  loss  mechanism  in  the  F-layer  has  been  identified  by  Dalgarno 
and  Deggs  [1968].  This  energy  loss  mechanism  arises  from  transitions 
among  the  fine  structure  levels  of  atomic  oxygen  and  is  induced  by 
electron  impact.  We  assume  that  energy  loss  to  ions  proceeds  by 
classical  elastic  scattering.  Schematically,  we  can  represent  the 
loss  term  by 


2  N  vei6et<69""  68+)  +  2  N  veo6eo6*’+  2  N  <v6)atonlc66' 

oxygen 


(10) 


where  the  last  term  represents  the  energy  loss  due  to  atomic  oxygen  . 

The  quantity  v  is  the  electron-ion  collision  frequency,  v  is  the 
electron  neutral  collision  frequency,  6.,  ®nd  6  are  the  average, 
fractional  energy  loss  per  collision  for  electrons  colliding  with  ions 
and  neutrals,  respectively  end  60—  is  the  departure  of  the  ion  (electron) 
temperature  from  ambient. 

The  quantity  Q  denotes  the  energy  source  due  to  the  RF  wave  heating. 
Q  depends  on  the  flux  of  energy  in  the  wave, 


F  .  J5L 
r  2  ’ 
4ttX 


(ID 


where  P  is  the  power  of  the  transmitted  energy  and  G  is  the  antenna  gain, 


4  _  16  D2  16  B2C 2 

G  -  ~2  ~~2 - 5 

A  C 


(12) 


*See  Appendix  for  interpolation  formulae  representing  physical 
processes. 
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D  is  the  dish  diameter,  £  the  transmitter  frequency  and  cp  the  half 
angle  of  the  beam  width. 

We  will  consider  energy  deposition  by  the  ordinary  ray  of  the 
magneto-ionic  waves  that  propagate  in  a  plasma.  For  the  longitudinal 
mode  (propagation  along  the  magnetic  field)  Q  is  approximately  given 


by, 


Q 


F 

c 


2 

%  v 
(u)  +  u^) 


(-efM 

cm  sec 


(13) 


where  is  the  usual  piasnu  angular  frequency,  y  is  the  electron 
collision  frequency  (assumed  small  compared  to  the  wave  angular 
frequency  uj)  , 

2 

UJ 

C  -  ! - (14) 

u>(u)  +  u^) 


and 


eg 
m  c 


(15) 


where  B  is  the  magnetic  field  strength  and  the  attenuation  A  is  given 
by 


A  “  exp 


r- 


,  X  -j)  2  u  dX 

I  f  .  P  --  a 

c  J  TT”  75  * 

o  s  (uj  +  a^) 


(16) 
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It  is  evident  that  for  high  frequencies.  u>  »  tu^,  the  heating  rate 
is  essentially  independent  of  frequency.  However,  as  u?  approaches 
the  penetration  frequency,  an  enhancement  in  the  heating  rate  takes 
place . 

For  the  ion  gas  we  write 

W  +  P+  M  ^  ■  I  »  »„•.!<«'-  ««+> 

m  N 


|  H  vl0  60+. 


(17) 


This  system  of  partial  differential  equations  was  converted  to 
a  system  of  coupled  finite  difference  equations  using  the  standard 
techniques.  The  heat  conduction  was  treated  explicitly. 

A  series  of  exploratory  calculations  were  run  for  various 
radiowave  frequencies  and  for  a  frequency  which  was  continuously 
adjusted  to  be  at  one-percent  above  the  frequency  which  nulls  the 
dielectric  constant  «  as  given  by  Eq.  (14).  The  results  of  these 
calculations  are  presented  in  the  next  section. 
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3.  NUMERICAL  RESULTS 

In  the  numerical  results  to  be  described  we  assumed  a  neutral 

composition  as  given  by  the  C1RA  1965,  Model  5,  hour  12.  Figure  1 

shows  the  neutral  composition  and  temperature  as  a  function  of  altitude. 

The  electron  density  as  a  function  of  altitude  is  ~hown  in  Fig.  2 

[W,.  H.  Hooke,  ESSA,  Private  Communication].  In  the  calculations  we 

assume  that  a  megawatt  of  average  power  arrives  at  an  altitude  of 

150  km.  Under  daytime  conditions  this  implies  that  considerably  more 

power  is  launched  from  the  ground  due  to  D-region  absorption.  Under 

nighttime  conditions,  D-region  absorption  is  essentially  absent  and 

although  the  neutral  particle  densities  are  appropriate  to  noon 

conditions,  the  nighttime  run  of  the  concentrations  are  at  times  typical 

to  those  shown  in  Figs.  1  and  2.  Thus,  we  may  consider  the  calculations 

to  be  presented  as  representative  of  conditions  resulting  from  a  one 

megawatt  average  power  ground  transmitter  with  an  assumed  aperture  of 
U 

10  square  meters  under  nighttime  conditions.  For  simplicity  we  assume 
the  electron  and  ion  temperatures  are  equal  to  ambient,  an  assumption 
more  appropriate  to  nighttime  conditions. 

In  the  numerical  work  we  took  the  ionosphere  as  extending  from 
150  km  to  800  km.  For  boundary  conditions  we  assumed  ambient  temperature 
and  pressure  were  maintained  on  the  top  and  bottom  of  the  ionosphere. 
Ten'kilometer  zone  sizes  were  assumed. 

The  relaxation  times  for  electron  cooling  due  to  collisions  with 
ions,  neutral  molecules  and  atomic  oxygen  are  shown  in  Fig.  3  under 
various  assumptions.  In  a  series  of  exploratory  calculations  done 


for  the  purpose  of  investigating  the  radiowave  frequency  sensitivity, 

we  assumed  electron  energy  loss  by  elastic  collisions  to  atomic  oxygen 

and  assumed  a  constant  fractional  energy  lo3S  per  collision  with 
-3 

molecules  at  10  The  cooling  times  under  these  simplified  assumptions 
are  also  shown  on  the  figure.  The  electron-neutral  collision  frequency 
was  taken  from  Mentzoni  and  Row  [1963].  Also  shown  on  the  figure  ate 
the  cooling  times  for  electron  cooling  by  rotational  excitation  of  • 
molecular  nitrogen,  fine  structure  transitions  in  atomic  oxygen 
and  electron-ion  cooling  assuming  elastic  collisions  with  atomic  it  ns . 
The  Appendix  provides  a  list  of  interpolation  formulae  used  to  represent 
the  various  physical  processes. 

As  stated  previously,  a  series  of  exploratory  calculations 
investigating  the  radiowave  frequency  dependence  of  the  plasma  heating 
and  expansion  process  were  performed.  The  results  are  shown  in  Figs. 

4  and  5.  The  maximum,  heating  takes  place  near  the  F-layer  peak. 

Figure  4  shows  the  fractional  change  of  the  electron  temperature  as 
a  function  of  time  at  an  altitude  of  300  km  and  for  frequencies  of 
10,  7,  6.35  MH*  and  for  a  wave  frequency  that  stayed  tuned  to  within 
one  percent  of  critical.  That  is,  to  »  frequency  such  that  c  ■  *01. 
The  required  frequency  change  as  a  function  of  time  is  indicated  on 
the  figure.  Evidently  the  electron  temperature  rises  to  a  plateau  in 
a  rather  short  time  (-*  1  min).  On  the  other  hand,  the  rate  of  expansion 
of  th*  plasma  is  quite  slow  as  indicated  in  Fig.  5.  We  attribute  this 

wjn  these  calculations  we  assume  s  magnetic  field  line  dip  of 
68° ,  appropriate  to  Boulder,  Colorado.  However,  in  order  to  assure 
penetration  of  the  wave,  we  assume  the  beam  is  directed  vertically 
upward. 


V. 
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to  the  drag  on  the  ions  induced  by  the  neutral  particles.  According 
to  this  model,  changes  in  electron  density  take  place  over  a  time 
span  of  tens  of  minutes  as  contrasted  to  nlectron  temperature  changes 
which  take  place  over  a  time  span  of  tens  of  seconds. 

As  a  final  exploratory  calculation  we  used  the  loss  mechanisms 
from  Mentzoni  and  Row  [1963]  and  Dalgarno-Degges  [1968].  The 
model  was  idealized  to  the  extent  that  we  assumed  vertical  field  line 
geometry  so  that  the  magnetic  field  and  the  axis  of  the  antenna  beam 
were  orthogonal  to  the  layers  of  stratification  of  the  ionosphere. 
Again,  we  varied  the  wave  frequency  to  be  within  one  percent  of 
critical . 

Figure  6  shows  the  fractional  change  in  electron  temperature  as 
a  function  of  time  for  the  indicated  altitudes.  In  all  cases  the  ion 
temperature  did  not  deviate  appreciably  from  the  initial  ambient 
value.  Figure  7  shows  the  electron  density  changes  with  time.  Al¬ 
though  the  numerical  calculations  were  not  run  to  a  final  steady- 
state  it  is  evident  that  reductions  in  electron  densities  of  the  order 
of  10  to  15  percent  may  be  expected  but  that  steady-state  conditions 
would  not  be  achieved  for  times  of  thirty  minutes  to  one  hour. 

An  interesting  feature  of  the  calculations  which  might  be 
susceptible  to  experimental  verification  in  the  predicted  behavior 
of  the  velocity  of  the  expanding  plasma.  When  the  beam  is  first  turned 
on,  a  wank  shock  is  generated  in  the  F-layer  peak  and  the  shock 
propagates  up  (and  down)  through  the  ionosphere.  Figure  8  shows  the 
velocities  generated  at  the  indicated  altitudes.  We  find  very  small 
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downward  velocities  (~  1  m/sec)  at  250  km.  (The  neutral  drag  inhibits 
the  motion.)  But  above  the  F-layer  peak,  where  the  neutral  drag  is 
much  smaller,  the  upward  velocities  reach  levels  of  20  to  60  m/sec 
as  indicated.  There  may  be  some  possibility  of  detecting  these  motions 
by  Doppler  techniques. 


14 


4.  DISCUSSION 

The  results  of  these  model  calculations  show  that  changes  in 

the  normal  F-layer  electron  density  by  factors  of  a  few  are  probably 

4  2 

not  obtainable  at  the  power-aperture  level  (10  Htu  *  m  )  to  be  used 
in  the  ESSA  experiment.  Changes  of  this  order  of  magnitude  have 
been  suggested  in  the  literature.  However,  if  the  Dalgarno-Degges 
[1969]  atomic  oxygen  cooling  mechanism  is  effective  then  one  will  be 
forced  to  higher  power-aperture  products,  even  when  following  the 
critical  frequency  to  within  one  percent. 

Another  interesting  feature  of  the  calculations  is  the  rather 
slow  time  scale  with  which  even  the  modest  10  percent  reductions  in 
electron  density  are  obtained.  This  is  due,  of  course,  to  the  viscous 
drag  by  the  neutrals  on  the  expanding  plasma.  This  treatment  of  the 
plasma-neutral  interaction  is  only  approximate.  However,  for  some 
forms  of  molecular  forces  this  form  of  the  drag  term  is  exact 
[Holway  1965].  In  any  event,  it  would  seem  that  one  must  be  prepared 
to  heat  the  ionosphere  for  times  of  the  order  of  tens  of  minutes 
before  substantial  reductions  in  electron  density  can  take  place. 

There  must  be  one  final  caveat.  We  have  assumed  non-diapersive 
propagation,!.?.,  the  effects  of  refraction,  as  the  electron  density 
changes,  have  been  ignored.  The  possibility  of  self-focusing  for 
horicontal  field  geometry  has  been  discussed  [Litvak  1968],  again  in 
the  steady 'state  approximation .  In  our  case  of  near  vertical  field 
geometry,  refraction  effects  may  be  expected  due  to  the  shape  of  the 
antenna  beam.  The  electron  heating  and  hence  the  rate  of  plasma 
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expansion  along  the  field  lines  is  greatest  at  the  center  of  the  beam 
and  decreases  as  one  moves  away  from  the  beam  axis.  This  effect 
produces  a  lowering  of  the  electron  density  which  essentially  reflects 
the  antenna  beam  shape.  It  is  easily  seen  that  an  electron  density 
of  this  structure,  lowest  at  the  center  and  increasing  as  one  leaves 
the  beam  axis,  leads  to  a  self "focusing  of  the  electromagnetic  energy. 
Crude  estimates  suggest  that  this  could  he  an  important  effect.  The 
GSSA  experiment  will  provide  both  in  situ  measurements  of  ionospheric 
electron  cooling  rates  as  well  as  a  possible  demonstration  of  self- 
focusing  of  electromagnetic  waves  in  the  ionosphere. 


v_ 
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APPENDIX 

FORMULAE  FOR  PHYSICAL  PROCESS 


Electron-ion  collision  frequency: 

-5  3/2  -1 

v  *  3  10  N/e  "  sec  (0  is  the  electron  temperature  in 

electron  volts) . 


Electron-neutral  collision  frequency: 

v  *1.95  lo”7  N  0  sec"*  (N  is  the  neutral  particle  density). 


Ion-  neutral  collison  frequency: 

.o  4.  + 

V.  -  6.4  10  N  g  sec  (g  is  the  ion  temperature  in  electron 
io  °  volts). 


Fractional  energy  loss  for  electron-molecule  collisions: 

6  «  (g  is  the  electron  temperature  in  electron  volts. 

eo  g  This  formula  is  a  fit  to  the  data  in  Mentsoni  and 
Row  f 1 963 ] • 


Loss  rate  to  atomic  oxygen: 

■  *•*  10'u  c°v«. 

oxygen 

(g  is  the  anient  air  temperature  in  electron  volts  and  [01  is  the 
a  atomic  oxygen  concentration.  This  is  a  rough  fit  to  the  data  in 
Dalgarno  and  Degges  [19681. 
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FIGURE  CAPTIONS 


Neutral  composition  as  a  function  of  altitude. 

(CIRA  1965,  Model  5,  12  hours). 

Ambient  electron  density  as  a  function  of  altitude. 

Electron  cooling  times  as  a  function  of  altitude. 

Fractional  electron  temperature  change  at  300  km  altitude 
as  a  function  of  time  for  the  indicated  frequencies. 

Fractional  electron  density  change  at  300  km  altitude  as 
a  function  of  time  for  the  indicated  frequencies. 

Fractional  electron  temperature  change  as  &  function 
time  for  the  indicated  altitudes.  The  radiowcve 
frequency  was  continuously  kept  within  1%  of  the 
critical  frequency  as  indicated. 

Fractional  electron  density  change  as  a  function  of  time 
for  the  indicated  altitudes.  The  radiowave  frequency 
was  continuously  kept  within  IX  of  the  critical  frequency. 

Velocity  cf  the  expanding  plasma  as  a  function  of  time 

for  the  indicated  altitudes.  The  radiowave  frequency 

was  continuously  kept  within  IX  of  the  critical  frequency. 


ALTITUDE  (KM) 


Fig. 3— Electron  cooling  time  as  a  function  of  altitude 
for  the  indicated  processes 


TIME  (SEC) 


Fig. 8— Velocities  induced  in  the  ionosphere  due  to  heating  process 


